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Abstract 

Supersonic  nozzle  flows  of  a  condensible  vapor  are 
considered  in  the  high  activation  limit  homogeneous  nuclea- 
tion.  Conditions  under  which  the  final  collapse  of  the  super¬ 
saturated  state  is  described  by  a  condensation  shock  are 
determined.  The  analysis  leads  to  a  natural  definition  for 
the  shock  position.  it  is  shown  that  the  shock  zone  is  asso¬ 
ciated  with  droplet  growth.  Droplet  production  occurs  upstream 
of  this  region.  Some  useful  similarity  laws  are  obtained  for 
the  production  zone. 
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1 .  Introduction 

t* 

In  two  recent  papers  non-equilibrium  nozzle  flows  of  a 
condensible  vapor  were  considered  in  the  limit  when  the  collapse 
of  the  supersaturated  state  occurred  close  to  the  condensation 
point  .  This  paper  is  concerned  with  flows  for  which  the 

collapse,  characterized  by  some  point  x  ,  occurs  at  a  finite 

distance  downstream  of  x  .  The  discussion  is  restricted  to 

c 

homogeneous  nucleation. 

for  this  process  the  rate  equation  contains  two  basic: 
parameters,  a  droplet  growth  rate  X  and  an  activation  factor 
K  .  it  appears  fShih,  1972)  that  A>>1  over  a  broad  range  of 
reservoir  conditions.  The  limit  X-*00  was  discussed  in  1  and  II 
for  flows  in  which  x.  x  a  o(l)  .  This  latter  requirement  can 

K  C 

be  viewed  as  a  restriction  on  the  magnitude  of  K  ;  at  fixed 
values  of  A  ,  x,  -x  decreases  as  X  increases.  In  fact, 
available  data  indicates  that  K<<1  in  many  practical  situations 
(\Shih,  1972)  and  the  presen'  "aper  is  concerned  with  the  double 
limit  ,  K-+0  such  that  \“xc  “  0(1)  . 

Of  particular  importance  in  this  analysis  is  the  behavior 
of  the  activation  function  B(p  ,T)  when  pv  is  the  vapor 
pressure  and  T  the  temperature  (52).  The  minimum  value  of  B 
corresponds  to  maximum  droplet  production.  for  nozzle  flows, 
the  variation  of  B  is  usually  of  the  type  shown  in  figure  1.  . 
Immediately  downstream  of  the  condensation  poi •  -  the  initial  be¬ 
havior  of  8  is  governed  by  the  frozen  solution  with  the  conden¬ 
sate  mass  fraction  g^O  .  It  was  found  in  1  and  II  that  the 
subsequent  growth  of  g  leads  to  a  strong  coupling  between  the 
rate  and  flow  equations.  This  coupling  induces  a  local  minimum 
in  B  as  shown  in  figure  1.  The  collapse  point  corresponds 

closely  with  this  turning  point.  It  is  important  to  note  that 
the  frozen  function  b„  also  has  a  turning  point,  at  x=x  .  , 

X  inin 

which  is  independent  of  the  parameters  X  and  K  .  For  some  of 

^Petty,  Blythe  and  Shih  (197?) >  Blythe  and  Shih  (1973),  referred 
to  respectively  as  I  and  II. 

Preceding  page  blank 


the  analysis  it  is  necessary  to  distinguish  between  flows  for 

which  x.  <xw,  and  flows  for  which  x.  ~x  .  .  (It  is  assumed 

k  min  k  min  v 

throughout,  consistent  with  standard  results  for  B(pv,T)  ,  that 

x  .  -x  »  0(1)  .) 
min  c  v  J  * 

Initially,  g  is  exponentially  small  and  its  growth  can 
be  computed  using  the  frozen  solution.  This  calculation  leads 
to  a  simple  criterion  for  the  determination  of  x^  (S3).  Qual¬ 
itatively,  is  strictly  associated  with  droplet  production. 

Although  the  production  rate  has  a  local  maximum  in  this  region, 
the  mass  fraction  g  is  still  small  ( 0  (K) ) .  Within  this  zone 
the  equations  reduce  to  a  unique  similarity  form  (§4)  which  is 
identical  with  the  corresponding  result  obtained  in  II. 

Downstream  of  x.  droplet  production  is  exponentially 
small  and  the  process  is  dominated  by  droplet  growth  (§5).  For 
x.<x  .  it  is  found  that  the  growth  region,  which  lies  close  to 
x^  ,  is  governed  by  a  condensation  shock  of  finite  strength  (§6), 
in  contrast  with  the  limit  discussed  in  II  where  the  collapse  was 
governed  by  a  weak  shock.  The  width  of  the  shock  zone  is  0(K^3) 

A  discussion  of  the  shock  relations  is  given  in  56.  In 
standard  treatments  it  is  conventional  to  replace  the  full  equa¬ 
tion  of  state  by  the  perfect  gas  law  (see  e.g.  Wegener  and  Mack, 
1958).  This  leads  to  a  considerable  simplification  in  the 
analysis  since  the  shock  relations  are  then  equivalent  to  those 
governing  the  one-dimensional  flow  of  a  perfect  gas  with  heat 
addition  (Shapiro,  1958).  It  is  shown  in  §6  that  the  terms 
neglected  in  the  classical  approach  may  be  important  at  high 
Mach  numbers. 

For  these  flows  it  has  often  been  shown  that  there  is  a 
critical  heat  addition,  or  condensation  fraction,  beyond  which 
the  flow  will  choke.  The  determination  of  this  crtical  amount 
is  somewhat  more  complex  for  the  full  shock  relations  discussed 
here  (§6). 

An  additional  constraint  on  the  mass  fraction  uis«. 
defined  by  the  rate  equation.  It  is  shown  in  §55  and  6,  neg- 
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lecting  terms  0 ( K  )  ,  that  the  limiting  value  within  the  shock 
zone  is  governed  by  saturation  conditions.  Consequently,  the 
flow  will  choke  only  if  the  limiting  saturation  temperature  is 
greater  than  the  choking  temperature,  No  particular  attention 
appears  to  have  boon  given  in  the  literature  to  the  determination  of 
conditions  under  which  saturation  will  precede  choking.  A  precise 
description  of  this  boundary  requires  detailed  numerical  calcula¬ 
tions.  However,  it  ir,  shown  in  appendix  1  that  choking  will  not 
occur  if  the  flow  at  the  condensation  point  is  supersonic  (Mc>l)  . 

As  xk*>'xmin  the  previous  remarks  require  modification. 

The  droplet  production  zone  now  depends  on  a  single  parameter 
v  .  Numerical  solutions  of  this  equation  are  given  in  §8.  Pro¬ 
vided  that  K  is  not  too  small  this  zone  is  still  terminated  by 
a  condensation  shock,  though  the  shock  width  is  now  0(K  6)  . 

For  still  smaller  values  of  K  (§9)  droplet  production 
is  governed  completely  by  the  frozen  solution.  Although  the 
rate  equation  does  simplify,  the  growth  zone  is  not,  in  general, 
described  by  a  condensation  shock. 

It  has  been  pointed  out  (Wegener  and  Mack,  1988)  that  the 
classical  shock  analysis  does  not  predict  either  the  shock 
position  or  the  end  state  immediately  downstream  of  the  shock. 

Nor  does  it  provide  information  on  the  conditions  under  which 
the  collapse  of  the  supersaturated  state  is  adequately  described 
by  the  shock  relations.  In  the  present  limit,  the  analysis  not  only 
yields  a  criterion  for  the  shock  position  but  also  outlines 
the  conditions  under  which  such  shocks  occur.  As  noted  above, 
if  the  collapse  is  governed  by  a  condensation  shock,  the  down¬ 
stream  limit  of  the  shock  zone  corresponds  to  a  saturated  state 

for  M  >1  . 

o 
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2 .  The  conservation  laws  and  the  rate  equation 

For  the  convenience  of  the  reader  the  relevant  results  are 
summarized  below.  These  relations  were  reviewed  in  I ,§2  where 


the  formulation  given  in  Wegener  and  Mack  (1958)  was  closely 
followed.  Under  the  assumptions  discussed  in  I,  the  equations 
governing  the  conservation  of  mass,  momentum  and  energy  for  a 
mixture  of  a  condensible  vapor  and  an  inert  carrier  gas  are 


puA  =  m  (2.1) 

<>u37*-3l  t2'2) 

and 

c  T  +  i  u2  -  g  =  c  +  i  u2  =  c  T  .  (2.3) 

P0  1  P0  1  c  P0  ° 

All  variables  are  normalized  with  respect  to  conditions  al  the 

4* 

condensation  point.  p  is  the  density  of  the  mixture,  u  is 

the  gas  speed,  p  is  the  total  pressure,  T  is  the  temperature, 

g  is  a  weighted  mass  fraction  (I,  2.12),  c  is  the  specific 

P  o 

neat  at  constant  pressure  for  the  mixture,  m  is  the  mass  flow 
rate  and  A(x)  is  the  local  cross-sectional  area  at  any  station 
x  .  The  suffix  c  corresponds  to  conditions  at  the  condensation 
point  x=l  where  p  =  p  =  T  =  1  and  g=0  . 

If  pv  denotes  the  partial  pressure  of  the  vapor  and  p. 
the  partial  pressure  of  the  carrier  gas 


i1  -  “.)  pi 


+  —  a)  p 
u  o*v 


(2.4) 


where  pQ  and  uy  are  the  molecular  weights  of  the  mixture 
and  the  vapor  respectively,  u  is  the  reservoir  specific 
humidity . 


S  defined  by  the  intersection  of  the  in i t i a 1  con- 
:  nil"  the  co-existence  line  . 


In  the  present  notation,  the  equation  of  state  for  the 
mixture  is 


p  -  p(l"H*1gjT  (2.5) 

where  H  is  a  non-dimensional  latent  heat.  (This  factor  arises 
m  (2.S)  due  to  the  normalized  form  or  g  ,  see  (2.oj  and 

I,  (2.12).) 

The  rate  equation  for  homogeneous  nucleation  can  be 
written 


where 


g 


fXM(x,OJ(f:)A(C)dC 


J  1 


M  (x 


•«  •  a: 


F(sj  (TD(s)-T(s)]ds 


(2.6) 


(2.7) 


is  the  mass  at  x  of  a  droplet  formed  at  £  and  J(x)  is  the 
rate  of  formation  of  droplets  of  critical  size.  It  is  assumed 
here  that  the  growth  function  F  is  independent  of  droplet  size 
with 


F  -  Ffpv,T,g)  .  (2.8) 

At  the  condensation  point  F  »  F(1,1,0)  =  1  .  In  (2.7)  T  is 
the  droplet  temperature  which  again  is  normalized  such  that 
TDal  at  the  condensation  point.  1’D  3  TD(p  ,T,g)  is  also 
assumed  to  be  independent  of  droplet  size.  It  is  important  to 
note  that  on  any  equilibrium  path  T=TD  » 

TD  =  Ts 

where,  in  the  size  independent  limit,  the  saturation  temperature 
is  defined  in  terms  of  the  local  vapor-pressure  by  the  in 
verse  of  the  Clausius -Clapeyron  law,  i.e. 

T"1  »  1  +  H ’ 1 1 n ( p  1 )  .  (2.10) 

s  v 
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The  droplet  production  rate 


J  -  £  (pv,T)exp{  -K~ 1  B(pv,T) }  .  +  (2.11) 

Here  B  is  proportional  to  an  activation  energy.  Classical 
theories  predict  that 

B  -  8  C  P  v  *  T ) (In  s)'2  +  (2.12) 

where 


s 


p  /p 
‘V  f! 


(2.13) 


is  the  saturation  ratio.  The  saturation  pressure  p  is  de¬ 
fined  by  the  Claus ius -Clapeyron  low 

pe  -  exp{  -H(y-- 1 )  }  .  (2.14) 

In  (2.12)  it  is  assumed  that  the  term  (In  s)  '  completely 
describes  the  singular  behavior  near  s-1  .  The  variables  are 
normalized  such  that  8(1, 1  “  1  . 


All  standard  models  «* 

If  necessary,  th**  analysis 


g  .  Voider,  1  9  V>  )  n  r  i-  f  t  h  1 1. 
can  be  generalized  to  i..  i  ..  .i  *- 


dependence  of  8  and  Y.  on  g  (sec  I  ,  l «  .)  . 


f  l  .  r  Ci . 

*4 


.T— 'it* 


^w^r.™t-!-nflr™wr  ^ww™r 


■j?  f*v^wjf?a 


* 


3 .  Initial  growth 

Since  B  is  unbounded  at  the  condensation  point  it  follows 
that  initially  g  !s  exponentially  small.  Tlu*  so'utlon  is 
determined  from  cho  frozen  relations 


-I_ 
T  f-1 

1  is 


*  0r 


JL.. 

T 
1  f 


and  u^,  •* 


rT  (W 


(3.1) 


whe  re 


c  T 

Po 


i  3 . 2  j 


This  system  is  completed  by  the  continuity  equation  (2,1).  IT 
the  condensation  point  li;„  downstream  of  the  throat  the  mass 
flow  is  defined  by  its  frozen  value 


m. 


/ 2T 


+  1 
H  MMI  I 

y-i 


A. 


(3.3) 


where  At  is  the  throat  area.  (3.3)  is  also  a  valid  firs*, 
approximation  when  no  appreciable  condensation  occurs  upstream 
of  the  throat.  If,  however,  the  supersaturated  state  collapses 
for  x<xt  the  mass  flow  is  determined  from  local  conditions  at 
the  appropriate  sonic  point. 

From  (2.6)  it  appears  that  the  initial  growth  of  the  con¬ 
densate  mass  fraction  is  defined  by 


A’f  r*  (x.OMOA^Jexpf-K  B.UJlde 
1  1 

•  [i  +  °(b/K)J 


(3.4) 


where  the  suffix  f  denotes  evaluation  from  the  frozen  solution. 
As  K-+0  this  rogion  of  near-frozen  flow  can  apparently  exist 
for  a  finite  distance  downstream  of  the  condensation  point, 
liquation  (3.4)  is  the  appropriate  extension  of  (II.  3.1)  and 
holds  even  when  the  initial  growth  region  is  not  confined  to 
the  neighborhood  of  x*=l  . 

11 


a  ..tv  ivr*ung  r-  |»  cv  usf^:  w»T,f*v*T*rv*'’r  u 


Near  the  eondonsat  ion  iio t t\t. 


T"l»  X  *«  •  S< 


(x-n 


U>  M 


apart  from  some  constant  factor.  This  limit  was  discussed  in 
li.  If  the  frosen  solution  persists  it  is  not  difficult  to 
show  from  dUftsic.nl  the  'ey  (1.4.8)  that 


H 


f 


l  .'.6) 


as  rf-0  ,  where  again  a  constant  factor  is  not  included  in 
(3.6).  The  characteristic  behavior  of  this  function  is  shown 
in  figure  l.  Only  situations  for  whicn  It  has  a  single 
turning  point,  at  say  x«xqUa  ,  will  be  considered  here. 

For  x>xmia  *  K+Q  ,  the  integral  in  (3,4)  is  obviously 
of  steepest  descents  type,  however,  the  cumulative  effects  of 
the  coupling  between  the  rate  and  flow  equations  can  become 
important,  upstream  of  *  In  this  case  the  coupling  can 

limit  the  droplet  production  and  lea. I  to  a  local  turning  point, 
character i ted  by  some  point  x»xk  »  for  the  activation  term 
8(PV,T)  (see  fig.  1).  The  analysis  givon  in  l l  was  concerned 
with  flows  for  which  *k-l  .  In  this  paper  the  analysis  is 
extended  to  situations  for  which  the  coupled  turning  point  can 
occur  at  a  finite  distance  downstream  -  f  the  condensation  point. 
It  is  convenient  to  distinguish  between  those  flows  for  which 
this  point  lies  upstream  of  xnl  and  those  for  which  it  lies 
close  to  xmltj  .  Initially,  the  analysis  is  restricted  to  the 
former  case.  Flows  in  which  ^ n  are  discussed  in  §  §  T - 9 . 

For  x<*oln  the  substitution 


t  ■  x  •  Ks 

can  ho  used  to  show  from  (3.4),  neglecting  the  factor  0(g/K)  , 
that 


J? 


g  •  1  Vl)}un'fU)A(K)i>*pt‘K'iBfU])  f 

’  » I 

•  II  ♦  om|  U,'*) 

provided  that 

K  ’  **  U;(kH  '  1  ,  (J,«) 

In  (3,7) 

...I'hLv  , 

rtj(x)  •  "  I  Pf  (OCT,*  fxM'fU)))  *  <*■») 

It  follow*  from  (A*?)  that 

8  •  bl  lKl|Uf  (x)E f  (x)  (8  |,( x) )  *  *V\p{  •  K  '  (» C -X ) ) 

•  f  l  •*  0  ( K )  |  .  (5,10) 

(5.4)  implies  that  this  initial  growth  law  will  tail  when 

g  “  O(K)  .  (5,11) 


lor  sufficiently  large  values  of  \  it  appears  from  i3.ld) 
that  (3,11)  will  occur  upstream  of  Aw,n  .  In  thin  case  it  is 
convenient  to  characterise  this  region  hy  defining  a  point  xK 
such  that 

— jJ-  .  D(xK)  ■  6A  '’k1j\*VV  (B£)  11  •  exp(  K  '  )  -  i  13.12) 

where  il  *  Q  ( x  )  etc.  As  implied  by  the  notation,  this  point. 

1  n 

will,  in  fact,  correspond  to  the  local  turning  point  discussed 
earl ier . 

The  criterion  (3.12)  may  he  of  some  practical  interest 
since  it  defines  a  critical  onset  point  for  the  collapse  of  the 
supersaturated  state.  Onset  criteria  huvo  been  discussed 
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extenalveiy  in  Ue  literature  (Wegener  (1900),  Wegener  and 
Mask  (1081)),  it  in  of  interest  to  observe  that  Onwatitaeh 
0^40  t  by  mean*  of  qualitative  a*«ueenta,  obtained  a  criterion 
which  in  similar  in  form  to  (3,420  In  addition  to  the  exponent  ini 
factor,  Oswatitach's  result  included  the  term  <  A 

similar  term,  (g--f)",l  arises  in  (3,12),  '* 

Strictly  the  point;  x^  defined  by  (3,12)  la  more  appro* 
pr lately  referred  to  at  a  relative  onset  point  for  which  the 
critical  value  of  g  ia  measured  with  respect  to  the  parameter 
K  .  Conventional  definitions  for  this  point  usually  correspond 
to  some  absolute  value  of  the  condensation  fraction  or,  at  leant, 
to  some  practical  measure  of  the  departure  from  «  frosr'n  state, 

For  noxtle  flows  th ’•  departure  is  often  defined  In  terms  of 
static  pressure  measurements  (see  e.g.  Wegener  and  Mack,  1958), 
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4 .  The  precursor  »one 

tournuiM#  M»wn»i  w*<iv  i -«  «r-ent-* 

Appropriate  asymptotic  axpanaiens  in  the  neighborhood  of 


XK  *»'!? 


etc . ,  where 


X  -  Kg, (♦j)  *  , . . 

T  *  Tk  ♦  ♦  ... 


(4.1) 


K  *(x-xk)  . 


14.2) 


these  expansions  are  compatible  with  the  limiting  behavior, 
x-*xk  ,  of  the  initial  growth  law  (3.10)  aiul  general  Ue  the 
approach  outlined  in  It, §4, 

Substitution  Into  the  conservation  equations  1,2 . 1 )  •  (2 . 3) , 
together  with  the  equation  of  state  (2.5),  gives 


O-O^n*.  *  (v-D(<-v',-'ikX4‘  . 


K-1'^  ■  •  *<Au».  *  t'  • 

•  \k*,  •  , 
K-Of’  •  -  . 


(4.1) 


Similarly,  it  is  not  difficult  to  show  that 

}  Pv  ,  I  2  -lT  2  (1-<0  Uv  2  VI  8, 

(Mk-D-‘  -  -  Y<*lk*,  *  |Ct-l)<-Hk,[YHj.  ^  Jj<  • 


In  (4.3)  and  (4.4) 


(4.4) 


(4.5) 
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«nd 


"k  *  "t;‘  • 


(4.6) 


Correspondingly,  the  rate  equation  (2.6),  using  (2.7)  and 
(2.11) ,  reduces  to 


r» 

r  (♦l“<>1),«xp{-B|1(i|»j))di|»  (4.7) 


81  " 


where 


R  -  \  ♦  KB^**)  *  ... 


and 


•-■(Si- 


(4.8) 


(4.9) 


From  (4.3)  and  (4.4)  it  follows  that 


with 


and 


whore 


B  j  «  -  ♦  bg 4  (4* 


a 


a.  •>('•*,).] 


b  '  N  (?  *  \  (P* 


\  ’ 


Tk(M*-!) 


(4.10) 


(4.11) 


(4.12) 


(4.13) 


::**n  ••  r?*r  T7y?ff  ? 
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In  this  section  it  is  assumed  that  xk<xratn  w,1^c^»  s*nce 
hj,  has  a  minimum  at  xmi  ,  implies  that  (tlBf/dx) ^*'0  or, 
from  (4*9) 


n  >  0  . 


(4.14) 


further,  for  the  usual  modols  (seo  e.g.  Vo  liner  1939) 


<  0  . 


(4.15) 


Only  flows  which  are  supersonic1  at  xk  will  be  considered  in 
this  paper.  From  (4.11)  to  (4.15)  it  is  then  not  difficult  to 


show  that  (M  >1) 


b  >  0  . 


(4.16) 


The  substitution 


$  »  a  4>  j 


(4.17) 


reduces  (4.7)  to 


g(<t>;b)  ■  [  (<p  -  ip)  3  exp  ( dJ  -  g  C  41  *.  b )  }  dip  .  (4.18) 


Although  the  mass  fraction  is  still  small,  0(K)  ,  (4.18)  implies 
that  the  coupling  between  the  conservation  and  rate  equations 
is  now  important.  The  region  is  characterized  by  a  marked  in¬ 
crease  in  the  droplet  production  factor  and  is  equivalent  to 
the  precursor  zone  discussed  in  II.  In  II  this  region  was  close 
to  the  condensation  point  but  here  it  occurs  at  a  finite  distance 


downstream  of  x  (=1)  . 

C 


The  solution  of  (4.18)  depends  on  the  parameter  b  .  How¬ 
ever,  the  further  transformation 


rWith  respect  to  the  frozen  sound  speed. 
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^v*~™**r 


yields 


♦  •  ♦  -  In  b 

A 

f4> 


g($)  "  i[  ($-$0  Jexp{$-g($)  )d$ 

u  J  *co 


(4.19) 


(4.20) 


which  is  parameter  freo.  (4.20)  is  identical  with  the  corre¬ 
sponding  equation  obtained  in  11,54.  Details  of  the  solution 
were  given  there.  As  $■*■»  it  was  shown  that 


g  -  as$*  +  aa$l  +  aj  +  afl  o(l) 


(4.21) 


where  the  error  term  is*  in  fact,  exponentially  small.  The 
coefficients  a,,  are  given  in  table  1  and  arc  defined  by 


/** 

5  J  3^r(**)3“P«xPf*-« (*))d*  .  (4.22) 


•T  » 

a 

0 

0.1883 

0  2220  1.0*73 

0.8560 

Table  1.  Tho  coeff‘ 


/fit*  o. 


From  (4.2)  the  width  of  the  precursor  tone  lit  apparently  0(K)  . 
It  is  relevant  to  note  from  ( 4 .  17)  that  an  effective  measure 
of  the  width  is  Kd  1  . 

Droplet  production  is  exponentially  small  downstream  of 
this  region.  The  cubic  growth  predicted  by  (4.21)  certainly  can 
not  persist:  this  growth  will  ultimately  be  limited  by  a  return 
towards  a  saturated  state.  For  tho  near  equilibirum  limit 
discussed  in  II  it  was  shown  that  the  final  growth  region  was 
governed  by  a  weak  condensation  shock.  The  corresponding 
analysis  for  the  present  limit  is  given  below. 
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Before  discussing  the  detailed  structure  of  the  growth 
region  it  is  convenient  to  re-write  the  basic  rate  equation  (2.6) 
in  an  alternative  form,  Using  (2.11)  and  the  definition  (3.12), 
the  integral  equation  (2.6)  becomes 

g(x;K)  -  exp{-K“1(B(f.;K)-B  )}dC  (5.1) 

6K  -'i  Hk  K 

where  the  implicit  dependence  on  the  parameter  K  is  now  in¬ 
cluded.  In  terms  of  a  normalized  radius 


R(x;K)  *  p  ds  , 


with 


x,  k 
k 


(5.2) 


SHx;K) 


F(Td-T) 


(see  (2.7)),  (5.1)  can  be  further  written  as 


g(x;K) 


~T  I  (-Dr  ,C  tR(x;K)]3'r  Ir(x;K) 

6K  raQ  3  r  r 


where 


(5.3) 


(5.4) 


a 

Ir(x;K)  -■  [R(g;K)]r  2I|.»KI  expI-K’1  (B(g;K)-Bk))dg  .(5.5) 

J 1  k  k 

Note  that  the  growth  of  R  is  measured  from 


x,  =  x.  -  Ka  In  b 
k  k 


(5.6) 


which  corresponds  to  the  displacement  defined  by  the  trans¬ 
formation  (4.19).  The  use  of  xk  ,  rather  than  xk  ,  enables 
the  precursor  solution  to  be  applied  directly  to  the  evaluation 
of  the  integrals  Ir  . 

Outside  of  the  precursor  zone  the  integrand  in  (5.5)  is 
exponentially  small.  From  the  precursor  solution  (§4)  it  is 
straightforward  to  show  that  for  x>x 
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Ir(x;K)  -  Ir(«;K) 


neglecting  exponentially  small  terms.  Further 


rxfl‘1’>r+1  r|,“> 

Ir(»;K)  -  l£-a^J- 

Hence,  from  (5.4),  x>x. 


exp[$«g($)]d$ 


g(x;K) 


a 


rF r  I  (Ka-l)r+1A  _  [RCxjK)]3^ 
bK  r*o  J 


where 


A,  _  -  (Ka‘1)'(r+1)  I„(-;K) 


3~r 


(5.7) 


(1+0(K) ]  .  (5.8) 


(5.9) 


(5.10) 


Using  (5.8)  and  (4.22)  it  follows  that 

Ar  -  &r  [l+0  (K)  ]  .  (5,11) 

Direct  inspection  of  the  equations  suggests  that  the 
precursor  expansion  fails  when  g«0(l)  or,  from  §4,  when 
x-xfe  *  O(K^)  .  Consequently,  appropriate  independent  and 
dependent  variables  in  the  growth  region  are  defined  by 


-  ZA 

x  -  xk  »  K  3  x 


(5.12) 


R(x;K) 


K^3  ft(X;K) 


with 


g(x;K)  -  g(X;K)  ,  T(x;K)  -  T(X;K) 


(5.13) 


etc.  Substitution  in  (5.9)  gives 

g(x;K)  -  £[a9a3ft3(x;K)  +  a2a2K^ft(X;K)  +  a^K^ftUjK)  +  0(K)] 


(5.14) 


where 


■ft(XJK) 


fXnrs ; 

*  n 


K)  ds 


(5.15) 


fl(X5K)  -  Qk1fi(x;K) 


(5.16) 


Obviously,  from  (5.15),  the  integral  form  of  the  rate 
equation  is  now  equivalent  to  the  first  order  differential 
relation 


fl(XJK) 


(5.17) 


with 


R(0;K)  *  0  . 


(5.18) 


The  function  fi  does,  of  course,  depend  implicitly  on  g  and 
hence,  from  (5.14),  on  R  .  This  relationship  is  defined  by 
the  local  solution  of  the  conservation  relations  together  with 
the  equation  of  state. 

Within  the  present  growth  region  the  dependent  variables 
have  asymptotic  expansions  of  the  form 


g(x;K)  -  g„ (X)  +  K  3gt  (x)  + 


(5.19) 


etc.  However,  it  follows  from  (5.12)  that  the  effect  of  the 

V 

local  area  variation  is  0(K  3)  and  that  the  conservation  relations 

reduce  to  the  standard  one-dimensional  form  even  if  terms  0(KXj) 

are  included  (see  below,  §6).  Correspondingly,  it  is  appropriate 

1^/* 

to  retain  terms  0(K  3)  in  (5.14). 

It  is  convenient  to  introduce  the  transformation 


s(x;K) 


'aa33[R(x;K)  +  ya"  K  3) 


(5.20) 


where 
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,-V,  V, 

X  ■  b  aa35  x  . 


(5.14)  now  becomes 

g(x;K)  *  ss(x;K)  ♦  o(KVs) 

with 


g(x;K)  -  g(x;K)  etc. 


In  addition,  from  (5.17) 


dS 

dx 


n(x;K) 


and  from  (5.18),  S  satisfies  the  initial  condition 


S(0;K)  -  a~2/i  (Kb‘l)  ^  ■  0.2253(Kb_1) 


■  K  Vj 


V. 


correct  to  0(K  3) 
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(S.21) 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


6 .  Condensation  shocks 

Substitution  of  the  expansion 

T(x;K)  »  T(x;K)  ,  g(x;K)  -  g(x;K)  etc.  (6.1) 


into  the  conservation  relations  (2.1)  to  (2.3)  shows  that  within 
the  growth  zone,  correct  to  0(K  3)  , 


pu  = 


p  +  pu' 


Cp  T  +  \  u2  -  i 
ro 


pkUk 


Pk  +  pkUk 


} (6-2) 


c  T.  +  \  u* 

p  k  2  k 


c  T 
P  o 
o 


where  the  right  hand  sides  follow  from  matching  with  the  pre¬ 
cursor  solution.  It  was  noted  in  §5  that  the  local  area  variation 
2/ 

is  0(K  3)  and  that  the  one  dimensional  relations  (6.2)  correspond 
to  the  usual  description  of  a  condensation  shock  (Wegener  and 
Mack,  1958). 


In  addition,  the  equation  of  state  (2.5)  becomes 


p  *  o(l  -H‘ *g)f  . 


(6.3) 


This  relation  is  often  replaced  by  the  approximation  (11-*^) 

p  =  pT  .  (6.4) 


Equations  (6.2)  and  (6.4)  are  equivalent  to  those  governing  the 
one-dinensional  flow  of  a  perfect  gas  with  heat  addition  (Shapiro, 
1958,  Wegener  and  Mack,  1958).  However,  neglect  of  the  term 
0(H  l)  in  (6.3)  can  lead  to  significant  errors  at  sufficiently 
high  Mach  numbers  (see  below).  It  will  be  retained  in  the 
present  discussion. 

The  solution  of  the  system  (6.2)  and  (6.3)  can  be  written 
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u 

°k 


p 


»1U0  ,YJ 

[y  +  1  *-(y-l)  gH’ 1  ] 

(6.S) 

where  the  stagnation  temperature 


(l  .  M’jl, 


1  ♦  Ml  M2 

Z  C 


(6,6) 


and  the  suffix  C  corresponds  to  conditions  at  the  condensation 
point.  Further 


£  , 
Pk 

1  ♦  YMi(l--  ) 
kv  u, 

k 

(6.7) 

and 

i  - 

rk 

■ 

(6.8! 

In  (6.5) 

the  function 

c  is  defined  by 

c(i 

H  A 

’V'J 

c'(vY)  f. 

■  c=(vY0Ji  -  ‘ 

:’(vY)(f.) 

3 

with 

(6.9) 

c-(vY) 

T 

*  v  ♦  i  -  2Y  g 

Y  1  t7!  rr  * 

(6.30) 

and 

Ci(vV) 

■  (2  *  Y  it')  n" 

iT  \  3 

(6.11 ) 

c.(vY) 

'  (Y'l)  (H°)  . 

( 6 .  i ) 

Consequently,  if  f  represents  any  of  the  dependent  variables 


u/u  etc.  then,  with  the  exception  of  the  vapor  pressure. 


f  *  , y)  .  (6,15) 

From  (2.4)  (I, (3. 2))  the  vapor  pressure  will  depend  on  the 


2 


.  1 

additional  parameter  P0rv  . 

Using  these  relations,  the  solution  of  (S.24)  can  be  written 

X  4  D  ■  f  rd-L  •  (*•**) 

• ^dO) 

where 

o(x;K)  -  ad(5)  16. is) 

neglecting  terms  GtK'**5)  .  The  dependence  on  the  parameters 
M  etc.  is  to  be  understood.  Using  (S.2S)  it  can  be  shown  that 

A. 

D  -  ^  a2a  j (Kb' 1 )  3  .  (0.16) 

Some  actual  calculations  of  the  shock  structure  for  a 
conventional  model  are  shown  in  figure  2  It  is  assumed  that 
Td*T3  (Oswatitsch,  1942),  w#»l  (pure  vapor),  and  that 

P  -  5*  p/  (6.17) 

(Wegener  and  Mack,  1958).  In  this  case,  including  the  parameter 
dependence, 

V  -  V(S;Mk,H,T0,y)  .  (6.18) 

For  the  calculations  shown  in  figure  2  the  downstream  limit  ol 
the  shock  solution  corresponds  to  a  saturated  state.  It  might, 
appear  from  the  rate  equation  that  this  asymptotic  limit  will 
always  be  defined  by  an  equilibrium  state  with 

T  ■  Tn  -  T  (6.19) 

I)  s 

(see  (2.9)).  However,  from  (6.5),  this  saturated  limit,  at 
which  g®gM  .  can  be  achieved  only  if 
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p  *  mJ  j*  c(^J?yv)  >  0  (6,ao) 


for  nil  gcgw  ,  It  should  be  atveased  that  g^  correaponda  to 
a  local  equilibrium  state;  it  is  not  defined  by  the  global 
solution  discussed  in  i,t,1. 

If  t>*Q  at  some  g*|w  the  flow  is  said  to  be  choked* 

It  is  straightforward  to  verify  from  the  shock  relations  that 
t,.a  flow  attains  sonic  speed  at  this  pcint,  where  the  appropriate 
sound  speed  is  defined  by 

*J  *  (Jp)#i<  *  nT>^fff7n  (5,a) 


and  the  su<  Uxes  denote  evsiuation  from  nn  Uentropic  frosen 
state. 

Although  it  is  important  to  understand  the  conditions  under 
which  choking  will  occur  upstream  of  the  equilibrium  state, 
surprisingly  little  analytical  work  ai peais  to  have  been  done 
on  this  problem,  In  moat  standard  treatments  the  classical 
shock  relations  (H*-**)  are  used  to  obtain  an  upper  hound  on 
the  mass  fraction  (heat  addition)  independent  of  the  equilibrium 
constraint.  It  is  shown  in  appendix  I  that  a  locul  saturated 
equilibrium  state  is  always  uttuined  upstream  of  choking  if 


X  ^  >  ix^JUL 

or,  equivalently, 


M  >  1  . 


Choking  is  possible  only  if 


M  <If 


or 


>  1  . 


(0,22) 


(0.23) 


(6.24) 


In  the  present  analyaia 


M  >  1 
k 


no  restriction  in  placed  on  Mc  . 
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£ 
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This  criterion  la  apparently  known  experimentally  (see  e«g.  Pour  inn  • 
1968 )  hut  it  does  not  appear  to  have  been  previously  deduced  by 
analytical  arguments, 

It  should  ha  stressed  that  (6,24)  is  a  necessity  rathav 
than  a  sufficient  criterion,  An  approximate  estimate  of  tha 
choking  barrier  is  deduced  in  appendix  i  under  tha  assumption 
that 

,U^1  (1«M*)  «  \  , 


it  it  shown  there  that  choking  win  occur  upstream  of  any 
saturated  state  if 


1>N‘ 


*  O  '  \ 

— >  C(H|Y»M»U.,  w. )  *  0  . 
(MJ’l)*  1  v  * 


whore 


2Cy-i)QCi^ Jff i ^ /i ♦  Q+ — ~ir~ 

mwwi—iw  —  ■*'iwmi'ir«  »*  m  wpig^MI  »A«phi  v>,  . . .  i  »  ivmi  wtwwSiA 


2(y4D  (Y*  U*Q* 


and 


(6 . 2S) 


(6,20) 


Q  *  1  yj y  *  pr  *  0  «  (0.27) 

More  detailed  numerical  calculations  are  presented  in 
figures  3-5,  !!or  given  H  and  y  these  figures  show  the  locus 
of  the  choking  and  equilibrium  points  in  the  (g,T)  plane  for 
a  pure  vapor.  The  numerical  results  confirm  that  choking  is 
possible  only  for  M  <1  and  (M  -1)  sufficiently  small. 

For  large  H  ,  the  shock  relations  simplify  considerably. 

In  particular 


C  -  (Y*l)  #  ,  (6.28) 

1  o 

and  the  approximate  choking  criterion  (6.25)  reduces  to 
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*  mo  1 

«MR«Min«  ^  ttJLv*»  ,  (^  t  29) 

(M*«l)*  yM 

However,  It  is  easily  seen  from  (0,0)- (6, It)  that  difficulties 
do  arise  with  the  Urge  W  approximation  when  T0  *  0(H) 

It  should  also  be  noted  that  for  the  finite  H  analysis 
given  here  the  choking  condition  0*0  way  have  more  than  one 
root,  In  fact*  if  choking  occurs ,  which  certainly  requires  that 
C^O*,  it  follows  that  the  cubic  (6.20)  has  two  positive  roots 
for  the  limiting  value  of  the  mass  fraction.  The  analysis  pre¬ 
sented  above  considers  only  the  smaller  of  these  two  roots,  in¬ 
dependently  of  the  saturation  constraint,  which  precludes  either 
root  for  M)s>i  ,  it  is  not  clear  whether  the  larger  root  corre¬ 
sponds  to  any  real  physical  situation, 


T 


Note  that 


for 


sufficiently  large 


C,  <0  . 
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7.  Initial  growth  |xk-xj  <<  J 

In  the  preceding  analyse  it  was  assumed  that 

»  1  (7,1) 

(see  (5.$)),  If  BJ»0  at  ,  (7.1)  can  be  replaced  by 

0  >  (V’Snin*  “  •  (7.2) 

u 

As  •  or  mor*  precisely  when  xit*xmin"°(K  )  *  an  alterna 

tive  approach  is  required.  The  frozen  turning  point  and  the 
actual  coupled  turning  point  are  now,  close  together  and  it  is 
convenient  to  Introduce  the  variable 

‘  •  •  <7-3> 

In  terms  of  c  the  initial  growth  low  (3.4),  neglecting 
the  error  term  0(g/K),  becomes 

g  ‘  >,K,nijAmoxp(-K'1Bffi)  [C  U-t)  Vxpl-|B"r*]dT 

*  m<Xi 

•  [1+OCK5)]  (7.4) 

where 

Bm  "  V*xmin>  etC*  (7*5) 

Corresponding  to  the  definition  (3.12)  it  is  appropriate  to 
define 


D  -  X  3KfJ*E  A  exp(-K'  *B  )  .  (7.6) 

m  m  m  m  © 

For  the  analysis  discussed  earlier,  with  8^0(1)  ,  Dm>>l  •  The 
criterion  (3.12)  implies  that  when  B]>-0(K1*)  ,  Dm«0(l)  .  Sub¬ 
stitution  of  (7.6)  in  (7.4)  gives 

g  -  KDm  [C  (C-x)5exp[-^T2]dr  .  (7.7) 

'  *00 


I 


V 
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It  is  apparent  that  (7.7)  is  not  valid  for  finite  c  when 
Dtt-0(I)  ;  the  orror  term  0(g/K)  is  then  important.  This  limit 
is  discussed  below  in  18. 

It  should  be  stressed  that  the  present  paper  is  concerned 
with  the  double  limit  X**.  ,  K+0  such  that  x^-1-0 (1)  .  At 
fixed  values  of  K  the  collapse  point  moves  downstream,  and  the 
parameter  Da  decreases  as  X  decreases.  The  analysis  for 
V°(i)  is  discussed  in  $9. 


..  www  i«rt --i 


8.  The  droplet  and  growth  zones,  Do"0(l) 


Tha  result  (7.7)  describes  only  the  initial  growth  of  the 
condensation  fraction  (c+-“)  .  Por  finite  values  of  c  the 
correct  expansion  is 


g(x;K)  -  g(c;K)  -  Kg2(c)  +  ... 

T  (x;K)  -  T(c;  K)  -  To  +  kVx(c)  +  KT2(?)  + 


}  (8.1) 


k 

etc.  It  is  easily  shown  that  the  terms  0 (K )  are  defined  by 
the  frozen  solution:  any  departure  from  the  frozen  state  is 
associated  with  the  terms  0(K)  .  Using  the  conservation  rela¬ 
tions  (§2),  it  can  be  shown  that 


(M2 - 1 ) T 
m 


l 


(Y - 1 ) (M*-y" r-H 


(8.2) 


with  similar  relations  for  the  remaining  dependent  variables 
(see  (4.3)).  Substitution  in  the  rate  equation  (2.6)  gives 


9a(C) 


a 


D 


m 


3exp [ - jB^t2 -bg  2 (t) ] dt  . 


The  transformation 


(8.3) 


9  -  bg2  ,  6  =  (jB^)15?  (8.4) 

reduces  (8.3)  to 
9 (c ;v) 

where 

v  -  4D  b (B") " 2  .  (8.6) 

m  m 


*  v  [  (5-r) 3exp[-T2-g(t;v)]dT  (8.5) 

J  -00 


(8.5)  is  the  appropriate  extension  of  (4.30)  and  governs  the 
droplet  production  zone.  The  simplicity  in  form  of  this  equation 
should  again  be  noted, though,  in  contrast  with  (4.10),  (8.5) 
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does  depend  on  a  single  parameter  v  . 

Numerical  solutions  of  (8.5)  were  obtained  for  various  v  . 

The  droplet  production  term 

3(£;v)  -  exp{-?2-fl(£;v)}  (8.7) 

i 

■  •  .  ■  .  'i 

is  shown  in  figure  6  and  the  normalized  droplet  number  density  I 

5  i 

W(c;v)  -  f  J (t ;v) dr  (8.8)  ’ 

is  given  in  figure  7.  In  addition,  the  asymptotic  level  N(«>;v)  , 

which  is  proportional  to  the  total  number  of  droplets  produced,  I 

is  displayed  in  figure  8.  Figure  9  shows  the  solution  for  the 

mass  fraction  g. 

The  asymptotic  growth,  ,  of  the  mass  fraction  is  again 

governed  by  a  cubic  law 

g  ~  oc3£3  +  ...  +  a0  +  o(l)  (8.9) 

where 

c*r(v)  3  3Cr  v  |  C-6)3"rexp[-e2-9(e;v)]dc  .  (8-.1C) 

These  coefficients  are  shown  in  figure  10.  Note,  in  particular 
that  a3(v)  «*  vN(«°;v)  . 

As  in  §5, it  follows  that  this  precursor  expansion  will  fail 
when  g“0(l)  .  The  asymptotic  law  (8.9)  implies  that  this  occurs 
when  c*0(K'^3)  or,  from  (7.3), 

*  ‘  *»i„  -  •  t8-U) 

It  is  easily  shown,  neglecting  terms  O(K^)  ,  that  the  dominant 
approximation  within  this  growth  zone  in  again  governed  by  the 
shock  relations  discussed  in  §6.  However,  the  error  term  due  to 
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the  area  variation  is^now  Q(K  *)  ,  significantly  larger  than 
the  contribution  0(1TS)  obtained  earlier  in  §§5  and  6  for 

VXmin  finite  C<0)  ’ 

Similar  results  to  those  outlined  in  §6  can  be  ^btnined 
for  the  shock  structure.  By  introducing  the  independent  variable 


b '  *  l/s  r  i.w *s v *  K 


'V  WK  '‘(x-xmln) 


(8.12) 


the  rate  equation,  within  the  growth  ione,  becomes, 


dS 

aw 


m 


neglecting  terms  OCK^)  .  Here 


(8.13) 


9  =  S 


(8.14) 


and  ftd($)  denotes  evaluation  from  the  shock  solution.  Matching 
with  the  precursor  solution  implies  that 


n 


rS 

O  Sd(4) 


(8.15) 
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i 

i 


] 

i 


9'  D*  :  o(i) 

It  was  noted  above  thac,  for  fixed"  K.  ,  the  parameter  D 
decreases  as  X  decreases.  Correspondingly,  the  collapse  point 
moves  further  downstream.  It  is  not  difficult  to  show  that  the 
behavior  near  is  described  by  (7.7)  with  an  error  0(KD2)  ; 

the  error  term  0(g/K)  in  (3.4)  does  not  play  a  significant 
role  for  •  o(l)  .  Droplet  production  is  governed  completely 
by  the  frozen  solution.  The  validity  of  (7.7)  is  actually  limited 
either  by  droplet  growth  or  by  the  local  area  increase. 

(7.7)  will  certainly  be  inappropriate  when  g  *  0(1)  ,  or 

C  -  0(KDm)‘^  (9.1) 

in  which  case 

*  •  Kin  ■  «(*V»d:S)  •  (9-2) 

Consequently  if 

KD;z  -  o(l)  (9.3) 

the  local  area  variation  will  not  be  important  within  this  region 
and  the  solution  is  again  governed  by  the  shock  relations.  In 

1/  _  Is 

this  case  the  shock  width  is  0(K  ®D  3)  and  the  error  term  is 
even  greater  than  that  obtained  in  §8  for  Dm=0(l)  . 

If  KDjj,*  is  not  small,  (7.7)  is  limited,  at  least  in  part, 
by  the  local  area  increase:  the  shock  relations  do  not  govern  the 
collapse  of  the  supersaturated  state.  Further,  it  does  not  follow 
that  the  flow  will  return  to  a  saturated  state.  Since  droplet 
production  effectively  terminates  upstream  of  che  growth  region, 
and  is  governed  by  the  frozen  solution,  it  can  be  shown  that  the 
rate  equation  will  reduce  to  the  simpler  form 
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_L 

with  3*0.  at  x"xmin  •  For  DffiK  *  0(1)  the  growth  law  (9.4) 
is  strongly  coupled  with  the  conservation  relations  and  the 
effect  of  the  area  variation  can  not  be  ignored. 


35 


Appendix  1 .  The  choking  condition 

From  the  supersonic  branch  of  the  shock  relations  it  is 
not  difficult  to  show  that 

—  >  0  (Al.l) 

ag 

upstream  of  choking.  The  initial  conditions  considered  here 
are  such  that 

Ts  =  (1+H” *ln  pv)'1  >  f  CA1. 2) 

and  >  0  .  It  follows  that  T  increases  monotonically 
through  the  shock.  Hence,  if  it  can  be  established  that 


f*  >  f* 

s 


(A1.3) 


where  *  variables  are  evaluated  at  the  choking  point,  at 
least  one  equilibrium  point  (T»T  «T  )  must  have  been  attained 

u  S 

upstream  of  choking. 


T* 


and 


The  shock  relations  imply  that 
2T6  (YM£M)2 


[i+Cy-1)|*] 


771  (Y*l)2M*[l  +  (^y)"(Mj~l)T 


jr~  (Al'.4) 


Y-l 


CM. 5) 


For  M*>1  it  is  straightforward  to  establish  that 
& 

Preceding  page  blank 
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(A1.6) 


and  similarly,  from  (A1.5) 


< 


(Y-l)f 


Obviously,  using  (A1.2),  if 


T£  >  [l+HMln(pS  r1]'1 

u 

the  inequality  (A1.3)  will  certainly  hold.  It  is 
that  (A1.8)  is  valid  if 


2T  o 
Y+T 


> 


1 


In  this  case 


ft >  1 

and  (A1.8)  follows  immediately  if  p*  <1  .  For  p 
more  care  is  required.  (A1.8)  is  equivalent  to 

exp[H(l-i  )  >  ?*  • 

'l  u 

Using  (A1.6)  and  (A1.7)  it  is  possible  to  re-write 


where 


and 


W,W2  >  1 


W. 


/  2T. 


(£t)  exp[H(l-^)_ 


w. 


[1-  Cy- 1 J 6 ] ' lexp 


(A1.7) 

(A1.8) 

shown  below 

(A1.9) 

(A1.10) 

*  >1  ,  a  little 
vu 

f  A1 . 11 ) 
(Al.ll)  as 

(A1 . 12) 

(A1.13) 

•'  A 1 , 1 4  ) 
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with 


(A1 .15) 


and 


fC«)  -  >  0  .  (A,  16) 

In  order  to  establish  (A1.12)  it  is  important  to  note  from 
the  shock  relations  that  choking  can  occur  only  if  Ct>0  or 


Y2- 1  >  22± 
y 


(A  1  . 1  7 ) 


If  (A 1.17)  does  not  hold  then  the  asymptotic  limit  of  the  shock 
solution  certainly  corresponds  to  a  saturated  state.  However, 
for  Cj>0  it  is  seen  that 


>  1 


wi  >  (ft)  ®xp[^t(f1  -1) 

provided  that  (A1.9)  is  valid.  Similarly 

.  1  r  Y  1  1  +  Y^T 

W2  >  [l+a-l)<5]  PXpJ^yffS)]  >  . 

Using  (A1.16)  it  is  now  straightforward  to  show  that 


(A  1 .18) 


( A 1 . 19) 


w2  >  1  . 


( A 1 .20) 


( A 1 . 18)  and  (A1.20)  establish  (A1.12)  subject  to  (A1.9).  Con¬ 
sequently,  the  asymptotic  limit  of  the  shock  solution  will  corre¬ 
spond  to  a  saturated  state  provided  that  T0  >  Ci;ii  o  r , 
equivalently  (see  6.6), 

M  >  1  .  ( A 1  .  2  1  ) 

c 

Conversely,  it  is  necessary,  though  not  sufficient,  that  Ml 
for  choking  to  occur. 
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Appendix  2.  An  approximate  choking  criterion 

The  condition  Mc<l  established  in  appendix  1  is  a  neces¬ 
sary  rather  than  a  sufficient  condition  for  choked  flow.  Numer¬ 
ical  calculations  (S6)  indicate  that  choking  takes  place  only 
if  Mk-1  is  sufficiently  small,  An  approximate  estimate  of  the 
actual  choking  barrier  is  given  here  in  the  limit  when  both 
Mr(>1)  end  Mp (< 1)  are  close  to  unity. 

It  is  convenient  to  introduce  small  parameters  c  and  A 
such  that 

Tq  ■  n**-^  (,1-c)  (A2.1) 

and 

M*  *  l  «  A  ,  (A2.2) 

Note  that  (A2.1)  is  equivalent  to 

c  •  ($x)(l-Mj)  .  f A 2.3) 

In  this  limit  the  amount  of  heat  required  to  choke  the  flow  is 
also  small  and  it  follows  from  t..he  shock  relations  that  at 
choking  g*  ■  0(A2)  .  Although  not  immediately  apparent,  it 
also  transpires  that  at  choking  c  “  0(A2)  . 

Expansion  of  the  shock  relations  shows  that  at  choking 

f*  »  1  -  e  +  — — ♦  [i».CY~p-Z|^  +  o (A 2 )  (A2..4) 

(7+1) *  l  Y+1  J 


and 

P*  -.1  -  “Xr  e 
*v  y  - 1 


yt*  ,  r yCv-13 
2Cy«1)!  L  y*1 


“/‘'Mi. 

"x".  J 


4* 


o  (  A  2  ) 


(A2.ii) 


where 


Preceding  page  blank 


6  =  g*/H  . 


(A2.6) 


further 


f»  « 

» 

i  *  ir‘((V* 

el*o  foiiaws  from 

i  <101  that 

IM  « 

<**.■**•*,  tJL V,  T,  f  1  , 

?(VT-1)  \ 

Choking  wiU  oceur  upstream 

'  \ 


i M  ?) 


(A3,H) 


f;  *  f* 

or,  using  tli#  preceding  relation*, 

V41  e  I’M* 

Y**  b  “  ;TiT7T)T  *  wf)  *  o 

where  B  is  defined  in  (o.aij). 


(A 2,S>) 
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Fig.  6.  The  droplet  production  rate  in  the  modified  precursor  zone 


Fig.  S.  Thf  Myaptotlc  4ropl»t  nu«b*r  u  a  function 
«(f  th«  pamtur  v  . 
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